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Visibility estimators and their performance are presented for use with the Palomar 
QO ' Testbed Interferometer (PTI). One operational mode of PTI is single-baseline visibility 

measurement using pathlength modulation with synchronous readout by a NICMOS-3 
H ■ infrared array. Visibility is estimated from the fringe quadratures, either incoherently, 

or using source phase referencing to provide a longer coherent integration time. The 
\q [ visibility estimators differ those used with photon-counting detectors in order to 

account for biases attributable to detector offsets and read noise. The performance of 
■ these estimators is affected not only by photon noise, but also by the detector read 

noise and errors in estimating the bias corrections, which affect the incoherent and 
coherent estimators differently Corrections for visibility loss in the coherent estimators 
O ■ using the measured tracking jitter are also presented. 
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1. Instrument configuration 



The Palomar Testbed Interferometer (PTI) flColavita et aL 1999| ; Wallace 1998 ) uses coherent 



fringe demodulation and active fringe tracking, similar to that employed with the Mark III 
Interferometer ( [5hao et al. 1988] ). Differences arise attributable to the use of an infrared array 
detector with its attendant read noise and required bias corrections. 

The beam combiner on PTI accepts the tilt-corrected, delayed beams from the two 
interferometer apertures. These are combined at a beamsplitter, and the two combined outputs 
directed to an infrared dewar. One output is imaged onto a single pixel of a NICMOS-3 infrared 
array. This white-light pixel is band-limited by an astronomical K (2.00-2.40 fim FWHM) 
filter, yielding an effective wavelength of ~2.2 [im. The other output is dispersed with a prism 
spectrometer and imaged adjacent to the white-light pixel onto the same line of the infrared 
array. Resolution is variable; one typical configuration uses 7 spectrometer pixels with center 
wavelengths of 1.993-2.385 /xm, yielding average channel widths of 65 nm. The combined light for 
the spectrometer channels is spatially filtered prior to dispersion with a single-mode infrared fiber. 
The white- light channel is not explicitly spatially filtered, although some filtering occurs because 
of the finite pixel size (40 /iin pixel and an f/10 relay). 
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2. Array readout 



The infrared array is read out coherently using a 4-bin algorithm with pathlength modulation 
implemented on the optical delay line. The 100-Hz modulation uses a sawtooth waveform, and the 
array readout timing varies according to the wavelength of each pixel to achieve a one-wavelength 
scan for all channels. Clocking constraints and overhead lead to a typical sample integration time 
of 6.75 ms (out of a sample spacing of 10 ms) for the white-light pixel, scaling proportionally for 
other wavelengths. 

For each sample period, the active and adjacent lines of the array are first cleared, the reset 
pedestal for each data pixel is read, and each pixel is then read out after each quarter-wave of 
modulation. Each of these (nondestructive) "reads" is actually an average of 16-64 consecutive 
2-jUs subreads, used to reduce the effective read noise, typically to a correlated-double-sample (cds) 
value of 12 e~ for the white-light pixel and 16 e~ for the spectrometer pixels. These 5 reads per 
sample for the white-light and spectrometer pixels are the fundamental interferometer data. 

Denote these 5 reads as Zi,ai,bi,Ci, and di, where i = denotes the white-light pixel and 
z = 1 ... 12 denote the R spectrometer pixels. The integrated flux in each quarter-wave time bin is 
calculated as Ai = a% — Zj, Bi = bi — ai, C% = c$ — bi, and Di = di — Ci. From these values, the raw 
fringe quadratures and total flux are calculated as 



The total flux Ni (as well as the integrated flux per bin) is related to the actual number of 
detected photoelectrons by Ni = krii, where A; is a dimensionless gain factor. For the PTI array 
electronics, k is typically 0.11. 

We can also calculate an energy measure which we denote as 



From these quantities we can estimate the fringe phase, visibility, and signal-to-noise ratio, but 
first it is necessary to correct for biases associated with the detection and readout process. 



Yi 
Ni 



Ai - Ci 
Bi-Di 

Ai + Bi + d + D { . 



(1) 
(2) 
(3) 



NUMj = X? + Y?. 



(4) 



3. 



Biases 



There are several bias terms that need to be measured. The first set of biases are the zero 
points of A, B,C, and D, and are those values observed with the instrument pointing at dark 
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sky.[](At high light levels, there are also nonlinearities as the detector saturates, but these effects 
are small for typical observations.) Expressed in terms of the quadratures and flux, we denote the 
biases as B X ,B Y , and B N , so that the corrected values of these quantities are given as (omitting 
subscripts for clarity) 

(5) 
(6) 
(7) 

We can also correct NUM for these biases as 

NUM* = NUM - B X (2X + B x ) - B Y (2Y + B Y ). (8) 

This is equivalent to simply computing NUM* as X 2 + Y 2 . 

The second set of biases occur in quadratic expressions like NUM and arise from the squaring 
of the photon and read noise. The two terms are just the variances 

B pn = kN (9) 

and 

B™ = Ak 2 a 2 cds . (10) 

Equation |9] is just the standard photon-counting bias. In Eq. a cds is the detector read noise 
(correlated-double-sample) , measured in the same units as the integrated flux. The factor of 4 
arises from the 4 bins used to compute NUM. We are usually read-noise limited on the channels 
of interest, in which case B rn dominates. Correcting NUM* for these two variances in addition to 
B X ,B Y , and B N yields 

NUM = NUM* - B pn - B Tn . (11) 



4. Bias measurements 

The biases for each pixel are measured at the beginning of each night of observation. While 
these initial values are adequate for proper operation of the real-time system, the biases are also 
measured repeatedly throughout the night for use in the science data processing. 



Initial calibrations A low-level calibration measurement is made at the beginning of each 
night with the instrument pointed at dark sky. The bias terms B X ,B Y , and B N are computed 
simply as the measured values of X, Y, and N. The bias term B rn is computed as the mean value 



1 With an ideal detector, these biases would be identical, proportional to the dark current and background. In 
practice, the biases on A, B, C, and D are slightly different. 
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of NUM*. This term also incorporates that fraction of the photon- noise bias attributable to finite 
dark count and background. 

A high-level calibration measurement is also made at the beginning of each night using 
an internal white-light source, which illuminates the white-light and spectrometer pixels. The 
increased value of NUM* with light level is used to estimate the pixel responsivity as 

k = (NUM* - B Tn ) /N, (12) 

so that B pn can be computed for other light levels using Eq. [9|. These values of B x , B Y , B N , B rn , 
and k are used by the real-time system. 

On-going calibrations Repeated measurements of the bias terms throughout the night 
accommodate drifts and improve the quality of the final data processing. Each typically 125-s scan 
on a science object is bracketed by several other calibration measurements: total-flux foreground 
and single-aperture ratio calibrations precede the scan; a background calibration, typically 25-s 
long, follows it. 

A foreground measurement observes the target with the instrumental pathlengths intentionally 
mismatched to yield zero fringe contrast. In this case, the observed value of NUM* can be used as 
a direct estimate of the sum B pn + B rn . The foreground calibration can also be used to estimate 
B x and B Y . 

A ratio calibration measurement observes the target with one aperture blocked. Combined 
with the total flux measured above, the intensity ratio between the interferometer arms can be 
estimated. 

A background measurement is essentially a low-level calibration measurement taken close in 
time to the stellar observation, and as such provides an estimate of B x , B Y , B N , and B Tn . 

These five calibration types can be used in different ways in the final data analysis. Typically, 
the biases B x , B Y , B N , and B in for each scan are estimated from the associated background 
measurement, while B pn is calculated from the actual flux during the scan using Eq. |9|. Averaging 
of several nearby background measurements using a median filter generally improves the calibration 
quality. The current data processing pipeline normally uses the foreground and ratio values only 
as diagnostics. 

5. Incoherent estimators 

Given the bias-corrected values X,Y,N, and NUM for the white-light and spectrometer 
channels, we can estimate fringe visibility. (Strictly, we estimate the square of the amplitude of 
the complex fringe visibility). Below we adopt a nomenclature for time intervals: a scan is a single 
measurement of an astronomical target, typically 120-150 s of recorded data, accompanied by 
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local calibration measurements as described above. A scan is divided into blocks, typically 25 s 
in length; the fluctuations of estimators among the blocks of a scan provides an estimate of their 
internal errors. Each block comprises a number of frames, which are typically 0.5 s long, and 
synchronized to the half-second tick. The significance of a frame is that intentional fringe hops to 
correct unwrapping errors in the real-time system are introduced only at frame boundaries. Each 
frame consists typically of up to 50 samples, which are data at the fastest rate in the system, 
typically 10 ms, which is of order of the atmospheric coherence time. The actual number of 
samples per frame will be less than 50 if fringe acquisition or loss occurs mid-frame; partial frames 
with less than typically 10 samples are discarded in the data processing. Squared visibility V 2 is 



estimated for each channel as (see Mozurkewich et al. 1991) 



v = ^2?2^. as) 

2 (AT)2 v ; 

where () represents an average over a block.0 While we're usually not photon-noise limited, the 
photon-noise-limited SNR is estimated similarly as 



The fringe phase is estimated as 



SNR2=2 252a. (14) 

(JV) 

(f) = tan" 1 — , (15) 

where we make no attempt to be rigorous with respect to phase offset. These estimates can be 
made for each channel: we typically use the suffix wl to refer to the white-light channel, viz. 
V 2 y = Vq. For the spectrometer channels 1 . . . R, we also compute a composite spectrometer 

7 2 

spec 

^spec ~T\ v~> / i\t \ 9tt7' ' (16) 



visibility V 2 C as 



it 2 E i <NUM i )W i 

^' 2 E*<^m 



The range of the summation covers channels 1 . . . R, or a subset (for example, 2 . . . (R — 1), which 
excludes the lower-flux channels at the band edges). The weights Wi can be uniform, but are 
typically computed as W{ = Nf/af^, wmcn are proportional to l/a 2 2 , as discussed below. 

3 i 

This composite estimator provides an improved signal-to-noise ratio, and is useful for compact 
sources where visibility changes with wavelength are smaller than the estimator noise. However, 
it still retains the wide fringe envelope (and thus decreased sensitivity to visibility errors caused 
by fringe-tracking errors) corresponding to the narrow spectral channels of the spectrometer; the 
use of the weights is useful for accommodating occasional spectrometer pixels with large read 
noises. For consistency, when the composite visibility is used for science, a composite wavelength 
computed with the same weighting is also employed. At the block level, the SNR of the V 2 
estimates is usually sufficiently high that the final V 2 estimate for the scan is calculated as a simple 
average of the block V 2 values, rather than carrying numerator and denominator separately. 



2 With step, rather than fringe-scanning modulation, the leading coefficient of Eqs. |l3| and ^ would be 4.0. 
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6. Coherent estimators 

We refer to the previous estimators as incoherent, in that NUM, the sum of the square of 
the fringe quadratures, is computed and summed over the 10-ms samples; these are generally our 
default estimators. However, the complex fringe visibility can be represented by a phasor; if the 
fringe phase is stable, we can add the phasors vectorially over multiple samples before computing 
NUM and related quantities. This "coadding" can provide an improved signal-to- noise ratio ( |Shao 



fc Colavita 1992a ), but at the expense of some atmospheric bias. To coadd the fringe phasors 
requires a phase reference, for which we use the white- light phase (p w \ = <pQ. 

We can compute a coherent visibility as follows: the white-light phase is scaled by 
the wavelength ratio between the white-light channel and the channel of interest to yield 
&i = 0wiA w i/Aj. The fringe quadratures are derotated and averaged over L samples as 

(^)coh = -^(Xicosfli-fisinfli) (17) 

(% oh = ^(Xisin^ + ficos^), (18) 

At PTI, the coadd time is typically one 0.5-s frame (L ~ 50 samples), although this is convenient 
rather than fundamental. A coherent value of NUM is computed as 

(NUM) coh = (X) c 2 oh + (Y) c 2 oh - (B pn + B™)/L, (19) 

where L reduces the bias correction to account for the reduced noise in the coadded quantities. 
From (NUM) co h and (iV) co h = (N) for each frame, the coherent V 2 is estimated as 

, = ^((N-UMU) 
2 <(iVU,}2 

A composite V 2 for the spectrometer channels can also be computed as for the incoherent case, 



similar to Eq. 16 



Given that the white-light channel has a high SNR, as required for real-time tracking, the 
coherent white-light V 2 is not an improved estimator because of coherence losses which occur 
in the phase-referencing process. However, it is valuable as an estimator of at least part of this 
coherence loss. We can estimate the coherence loss T a as 

It - (V^) coh /V 2 h (21) 

and we usually divide the coherent spectrometer V 2 values through by this value as a partial 
correction. To be more exact, one can account for the wavelength difference between the 
white-light and spectrometer channels by scaling the correction with wavelength as 

which assumes a simple exponential form for the coherence loss. We note that there are additional 



^-(fJMT- (22> 



coherence losses in phase referencing, some of which are discussed in Sec. 8.2 



-7- 



7. SNR of the V 2 estimators 

The "detection" noise on the V 2 estimator attributable to photon and read noise (as opposed 
to noise attributable to atmospheric turbulence) is readily calculated. As is usual, we model 
only noise on NUM, given by Eq. ^, and ignore the smaller noise in N that normalizes NUM 
in calculating V 2 ( Tango fc Twiss 1980| ). The quadratures X and Y are each comprised of two 



correlated-double-sample reads, so that the variances of X and Y are given as a\ = a\ = 2a 2 ds . 
For additive Gaussian noise, a 2 x2 = = 8cr 4 ds , and similarly for Y, yielding ctnum = 4<x 2 ds . 
Thus, the standard deviation of the (incoherent) V 2 estimate in the read-noise limit is 

2?r 2 /a cds x2 



av2 = 7W\J^ ] > N<<N ™> (23) 



where M is the total number of samples, both temporal and spectral, in the estimate, and is thus 
applicable to both single-channel and composite (with equal weights) visibility estimates.^ For 
arbitrary photon fluxes, read noise can be incorporated into the standard (4-bin) photon-counting 
result ( Tango fc Twiss 1980| ), yielding 



^ = 4^( Ar2 + > V2 + 16 ^)' <24> 
Thus the read- noise limit applies when N <C N m , where 

N rn = min (4a c 2 ds , 3AV~ 2 / 3 a^) . (25) 



A numerical example is illustrative. For the case of a read noise of 16 e~ per pixel, 125 s 
of data at 10 ms per sample, and 5 spectrometer channels in the spectral composite, a standard 
deviation of 0.02 requires 32 photons per channel per sample. 

For the coherent estimators, the standard deviation is similar. Assume as above that M is 
the total number of 10-ms samples in the estimate, but that they are first coadded to frames of 
length L before calculating NUM. In this case 

2vr 2 fa cds \ 2 



^ = 7m{-w) • JV<<JV - (26> 

Thus, the required photon flux for a given accuracy scales as L -1 / 4 . With L = 50 and the 
parameters above, an accuracy of 0.02 now requires 12 photons per channel per sample, although, 
as discussed above, the coherent estimate is more susceptible to systematic biases. 



3 With step modulation, the leading coefficient of Eqs. [27| and ^ would be 16.0, with similar changes to 

Eq. ^. 



- 8 - 



7.1. SNR for bias estimation 

Strictly speaking, the above analysis is somewhat simplistic, as it assumes that bias correction 
adds no additional noise. For low light levels, the largest errors in bias correction are attributable 
to estimation of the biases B pn and B rn ; errors in their estimation are additive with the noise on 
NUM as calculated above. However, as B pn and B rn are computed from NUM measured under 
known conditions (Sec. ||), the V 2 errors due to imperfect bias estimates can be computed using 
the expressions above. Thus, for incoherent quantities, the (incoherent) V 2 error due to imperfect 
bias subtraction is given by 

e v2 ^^=(^) 2 ,N«N ia , (27) 



HvT b V N 

where Mf, is the number of samples used in estimating the bias. This expression is strictly accurate 
only for the read-noise bias B rn , or when both B pn and B rn are computed from a foreground 
calibration. However at low photon fluxes, where bias errors are most significant, the read- noise 
term is dominant and the above expression is a good approximation. 

For the coherent V 2 , the situation is somewhat better, as the errors in B pn and B Tn are 
reduced by the number of samples in the coherent average, per Eq. For the biases computed 



incoherently, and applied via Eq. |19|, the applicable expression is 

,2 . \ 2 



subject to the same caveats at Eq. 27. By way of comparison, if the biases were computed 



"coherently", i.e., from measured values of {NU M) co ^, then Eq. 28 would have the same 



dependence on L as Eq. 26 



Thus, the total "detection" noise on V 2 is the quadrature sum of a and e, and the contribution 
due to bias estimation can be important. This contribution is generally not important on bright 
sources where the noise on V 2 is dominated by atmospheric effects. On fainter targets, the relative 
bias noise can be decreased by incorporating additional calibration data (for example, using 
background calibrations from a larger time window about the science scan, rather than just its 
explicitly-associated background), although eventual nonstationarity of the underlying statistics 
presents a practical limit. 



While calibration errors are usually dominated by the B pn and B rn terms, the errors 
attributable to the other bias terms are easily computed: for both incoherent and coherent 

X f Y anrl f N 

y2, ty 2 , ana ty 2 



estimators, the errors €y 2 , €y 2 , and Cy 2 associated with B x , B Y , and B N are given by 



* = * = <29> 
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8. Data quality measures 

Inter-block fluctuations of estimated quantities are useful to estimate internal errors. However 
additional data quality measures are available. 



8.1. Lock time 

PTI uses a multi-stage algorithm for fringe acquisition and track ( Polavita et al. 1999| ). 
Essentially, the average SNR must exceed a given threshold for the system to enter the "locked" 
state; loss of lock and reacquisition occurs if the SNR falls below a second threshold. Fringe data 
is only recorded when locked; to account for the time delay caused by the memory of the averaging 
filter in detecting loss of lock, data at the end of a lock is automatically expunged. Thus, with 
multiple locks, the elapsed time to collect a fixed amount of data in order to complete a scan is 
increased. 

Two heuristic data-quality measures are the fraction of lock time to elapsed time, and the 
number of separate locks that make up the total data on a scan. For bright stars and good seeing, 
each scan is comprised of just several long locks. For very faint stars, or with poor seeing, each 
scan can be comprised of many short locks, reflecting the inability of the system to consistently 
track the fringe. While visibility can be estimated in all cases, the data quality in the latter case 
will be inferior. Typically, this poorer data quality is evident in the inter-block fluctuations, in 
which case the lock-time metric is only advisory. 



8.2. Jitter 

We can estimate a first-difference phase jitter as 

°l<t>= <(<Mn)-<Mn-l)) 2 >, (31) 

where w i is computed from the 10-ms samples. While this quantity is not unbiased with respect 
to detection noise, successful fringe tracking typically requires an SNR > 5, so that the detection 
bias on should be < 0.08. 

With an ideal instrument, a&s is related to the atmospheric coherence time. Coherence time 



can be defined in various ways (Buscher 1994). Let io,i denote the structure- function definition of 
coherence time, viz. that sample spacing for which the phase difference between samples is one 
radian rms. The structure function depends on the actual sample spacing t as D{(t) = (i/to,i) 5 ^ 3 - 
For i = 1, representing the contribution from a single aperture (the usual adaptive-optics 
definition), io,i = OSIAtq/W for coherence diameter r$ and constant wind speed W; for % = 2, 
representing contributions from two apertures (applicable to interferometry) , io,2 = 0.207ro/W. 
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Let To t i to denote the variance definition of coherence time, viz. that sample integration 
time for which the phase fluctuations about the interval mean are one radian rms. It is given by 
T 0; i = 1.235r /W and T 0:2 = 0.815r /W. The structure function depends on the actual sample 
integration time T as of = (T/T QA ) 5 / 3 . 

Thus, for an ideal instrument, the coherence time Tq^ can be estimated as 

T 0i2 = 3.94t/a^ 5 , (32) 

where t is the sample spacing. Fringe motion during the sample integration time T blurs the 
fringe, reducing the visibility. For rapid (with respect to underlying phase motion) fringe scanning, 
the coherence reduction is related to the high-pass fluctuations about the interval mean, (cr^)hp, 
as T b = exp(— (<7^)hp), or given the coherence definitions above, T b = exp(— (T/To^) 5 ' 3 )- We can 
write this in terms of the phase jitter as 

T b = exp (-C r a 2 M ) , (33) 

with the coefficient Cr given by 

C ^(lMi) ' < 34 » 
For T = 6.75 ms (for the white-light pixel) and t = 10 ms, Cr = 0.053. 

A more careful calculation of Cr can be done for this case (Appendix |A|) , accounting for the 
finite integration time required to measure <p w \, which yields Cr = 0.057. A similar calculation 
can be done under the assumption that all phase noise is caused by narrow-band vibrations with 
frequency <C 1/t; in this case, Cr = 0.038. 

When we apply this correction to PTI data, we usually err on the side of undercorrection 
by adopting a modest leading coefficient of 0.04. In general, an empirical visibility-reduction 
coefficient can be adopted from fits to the measured data applicable to the actual atmospheric 
realization and instrumental configuration. However, for data calibrated with spatially- and 
temporally-local calibrators (and especially if the calibrators are of similar brightness to the 
target), the reduction in visibility due to the above temporal effects will be mostly common mode 
and divide out of the normalized visibility. In this case, the value of the jitter is useful as a 
measure of the seeing, and indirectly of the data quality. Finally, we note that the coherent V 2 
estimates on PTI often exhibit coherence losses larger than predicted from the models above. 
These may be attributable to different apodizations of the starlight pupil between the spectrometer 
and white-light sides of the beamsplitter. In particular, the single-mode fiber preceding the 
spectrometer imposes a Gaussian apodization on the pupil, while the white-light channel — with 
no explicit spatial filter — imposes a more uniform pupil weighting. These different apodizations 
will result in slightly different instantaneous phases between the two beamsplitter outputs, and 
thus a coherence loss when phase referencing the spectrometer channels to the white-light phase. 
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8.3. Ratio Correction 

PTI uses a single-mode fiber after beam combination to spatially filter the spectrometer 
channels. Spatial filtering increases the raw visibility and reduces the concomitant noise 
attributable to spatial effects; temporal effects must still be calibrated. As spatial filtering by the 
fiber essentially rejects light which would not interfere coherently, there is induced scintillation, 
which has a second-order effect on visibility. With simultaneous intensity measurements of each 



arm in a fully single- mode combiner ( Foresto 1994 ), an essentially perfect correction for this effect 



is possible, but it can be shown ( fghaklan, Colavita, fc Shao 199^ ) that measurement of only the 
average intensity ratio between the two arms is adequate. If we denote this ratio as R12, then the 
correction for the induced scintillation is 

(1 + R 12 ) 2 , . 

Sl2 = 4R 12 ■ (35) 

As discussed in Sec. [|, the combination of the foreground and ratio measurements allows estimation 
of S12 for each scan. 

Currently, strict application of the ratio correction at PTI has been unsatisfactory, and we 
generally do not apply it. We attribute this to two effects. One is that given noisy values of 
R12, S12 is a biased estimator, and will tend to over-correct the visibility. The second is that the 
measurements of the ratio are not truly simultaneous with the scan. Thus, seeing nonstationarity 
will affect the estimate. Also, there is a selection effect as fringe data is only recorded when locked, 
while the flux calibrations are contiguous. 

Even without the ratio correction, the spatially-filtered data yields significantly-improved 
visibility estimates. However, the ratio correction has been useful as an additional indicator of data 
quality. For example, at high zenith angles, asymmetric (due to misalignment) vignetting in the 
system will increase S12 ■ But as with jitter, vignetting is tied to sky position, and spatially-local 
calibration will ameliorate most of the systematic visibility effects. 



9. Conclusion 

The use of array detectors at PTI requires attention to bias correction in fringe-parameter 
estimators, especially energy measures like V 2 which use squared quantities. Observations with 
PTI incorporate nightly and on-going bias calibrations, which can be used to compute optimal 
bias corrections. In addition to statistical noise in the estimators themselves, noise in the bias 
terms plays a role in the overall data quality. Inter-block fluctuations of estimated quantities are 
useful to estimate internal errors. Auxiliary data quality metrics include the tracking jitter and 
the ratio-correction estimate, which can be used for open-loop corrections or as independent data 
quality measures. 
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A. Temporal coherence calibration using the phase jitter 

We start with the assumption that the coherence reduction on V 2 can be written as 

r = exp (-K)2 p ) , (Al) 

where (c</,)hp is the high-pass phase jitter. This is strictly true for the case where the fringe 
scanning is much faster than any frequencies of interest, although the results are similar with a 
slower scan. The high-pass jitter in Eq. Al is given by the frequency-domain integral 

a 2 p = J dfW(f)(l - sincV/T)), (A2) 

where W(f) is the phase power spectrum, 1 — sinc 2 () is a high-pass filter, and T is the sample 
integration time. A similar spectral representation exists for the phase jitter (Eq. |3l| ): 

°U = I dfW(f) S mc 2 (7rfT)4 S m 2 (7rft), (A3) 

where sinc 2 () accounts for averaging over the sample integration time, while sin 2 () is a high-pass 
filter corresponding to a sample spacing of t. 

For / < 1/T, the filter function in the integral for a^ p is H hp (f) ~ ^ir 2 T 2 f 2 , while for / < 1/t, 
the filter function in the integral for is i?A</>(/) — 4vr 2 i 2 / 2 . The ratio of the filter functions is 

C r .##.i(?V (A4) 



H A 4f) 12 V t 

With T = 6.75 ms (for the white-light pixel) and t = 10 ms, we calculate Cr = 0.038. Thus, for 
narrowband low-frequency noise, we can write the visibility reduction directly in terms of the 
first-difference variance as 

T b = exp (-C r ^) . (A5) 

This same formulation applies for other noise models. For W(f) given by an atmospheric 
power spectrum, W(f) oc /~ 8 / 3 (assuming a low fringe-tracker bandwidth), it is necessary to 
compute the integrals numerically. For power laws of the form f~ a , some representative values of 
C r for T/t = 0.675 are 0.057, 0.070, 0.088, and 0.145 for a = 8/3, 2.5, 7/3, and 2.0, respectively. 
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